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1.  Introduction 


As  a  result  of  the  December  1978  workshop  session  of  the  panel  on  gravity 
field  and  sea  level  requirements  the  resulting  report  "Applications  of  a  Dedicated 
Gravitational  Satellite  Mission"  (1979)  contains  a  discussion  on  the  kind  and  pre¬ 
cision  of  gravity  parameters  to  be  provided  through  a  dedicated  gravitational 
satellite  mission.  These  requirements  are  in  essence: 

1.  For  geological  or  geophysical  applications  gravity  data  (gravity  anomalies  ?) 
with  a  half-wavelength  resolution  of  down  to  100  km  and  an  accuracy  of  ±2.5 
to  10  mgals  depending  on  the  application; 

2.  Oceanographers  would  need  geoid  height  differences  for  distances  of  100  to 
3000  km  with  an  accuracy  of  ±10  cm. 

The  same  report  contains,  in  addition,  a  list  of  open  questions  concerning  e.  g. 
the  assumptions  on  the  range  rate  data  noise,  the  needed  orbital  accuracies,  the 
''best"  gravity  parameterization,  the  appropriate  error  measure  a.s.o. ,  to  be 
analyzed  in  further  studies. 

The  purpose  of  this  report  is  to  answer  part  of  these  questions  in  an  error 
analysis  for  satellifce-to-satellite  tracking  (SST)  in  the  low-low  mode.  For  the 
influence  of  systematic  and  random  uncertainties  of  the  satellite  orbit  upon  the 
gravity  estimates,  a  problem  of  high  relevance,  which  is  considered  in  Douglas 
et  al.  (1980)  too,  sensitivity  coefficients  are  derived.  The  study  will  also  analyze 
the  dependency  of  the  accuracy  of  the  estimated  gravity  parameters  on  the  meas¬ 
urement  precision,  the  altitude  of  the  satellite,  the  separation  and  spatial 
arrangement,  and  on  the  data  coverage.  As  a  conclusion  of  the  above  stated  user 
requirements  the  estimated  parameters  will  be  on  one  hand  l°x  1°  mean  gravity 
anomalies  (for  geophysical  purposes)  and  on  the  other  hand  geoid  heights,  and 
geoid  height  differences,  respectively,  (for  oceanographers).  The  estimation 
model  that  allows  without  difficulties  this  type  of  flexibility  is  the  least  squares 
collocation  model.  Because  of  some  conceptual  problems,  this  model  will  be 
simplified,  without  any  serious  impact  on  the  validity  of  the  results,  however. 
Special  attention  will  be  put  on  the  definition  of  a  valid  concept  for  the  error 
measure. 

This  study  may  be  considered  as  complementary  to  those  of  Schwarz  (1970) 
and  Kaula  et  al  (1978)  in  that  a  different  type  of  mathematical  model  is  employed, 
to  Krynski  (1978)  because  of  the  more  operational  type  of  model,  and  to  Rummel 
(1979)  where  fora  purely  theoretical  model  the  optimal  situation  is  analyzed.  It 
is  also  complementary  to  Rapp  and  Hajela  (1979)  who  discuss  In  a  very  similar 
manner  the  corresponding  SST  high-low  case. 
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2.  Observational  Model  for  Range  Rates 


The  range  rate  jb  between  two  satellites  at  Sj  and  is  expressed  as 

»Sa)  =  X is  *  »  (1) 

•  •  • 

where  Xls  =  X(S2)  -  X(Sj)  is  the  velocity  difference  between  the  two  space 
probes,  and  £.12  ~  0  X 13  •  the  unit  vector  pointing  from  (with  position  vector 
Xi  ,  expressed  e.  g.  in  an  quasi-inertial  frame)  to  S2  and  p  is  their  distance. 
From  equation  (1)  the  observational  model  will  be  derived,  similar  to  the  model  in 
Rummel  et  al.  (1978)  and  according  to  guidelines  developed  in  Grafarend  (1979). 
[Neither  the  original  character  of  the  range  rate  observations,  as  described  in 
Eddy  &  Sutermeister  (1975)  or  in  Barthel  et  al.  (1978),  which  may  be  for  example 
fringe  or  doppler  counts,  nor  the  effect  of  non-gravitational  and  tidal  perturbations 
on  the  motion  of  the  satellites  will  be  considered  in  this  context.] 

Each  observed  range  rate  at  time  ti  is  linked  to  the  gravity  field  of  the 
earth,  parameterized  by  £  =  £c  +  6£  (  3C . . .  known  computed  longwavelength 
part  of  the  gravity  field,  . ..  residual  part)  through  the  position  and  velocity 
vectors,  X  and  X,  of  the  two  satellites-. 

li  =  4(t1;  Xi(g)>  Xi(£).  X3(g))  +  e,  (2) 

with  fi ..  .observation  error.  With  Y:T  =  ( Xi ,  Xj. ) T  for  the  complete  state  vector 
and  with  Yl0  for  the  state  vector  at  t  =  to  (initial  condition)  a  Taylor  expansion  of 
equation  (2)  at  the  computed  reference  locations  S'  and  S2  of  the  two  satellites 
yields  j 


it  =  4*  (*i »  S' ,S2C)  +  (gradyj  f>)\  (S^ )  [(grady^  Y1)T(Sie)  AY10  + 

(gradg  YX)T  (S*)Ml  +  (gradya(5);  (SDUgradv^Y^)'  (S|)AYao  + 

(grad^ Y/  (Sae)6g]  +  f,,  (3) 

where  e,  shall  now  contain,  besides  the  measurement  noise,  the  terms  of  second 
and  higher  order  of  the  series  expansion.  [The  gradient  of  a  vector  shall  thereby 
express  the  gradient  of  the  components  the  vector  and  constitute  a  matrix  of 
dimension  (no.  of  gradient  components)  x  (no.  of  vector  components)].  Equation  (1) 
yields: 

fPt~1(X!a  - 

-  (gradIi  t»(S{)  =  (gradr,  (5)(S,C)  =  ^ 


and  in  analogy  to  Schwarz  (1970)  we  denote 


(gxadv^  Yi)  (S?)  =  ix , 
(&radIao  la)  <S«)  =  4a> 

(ii  .  I3lf  =  iT. 

(grad  g  Y,)  (Sf )  -  Wx  , 
(grad  Y,)  (S‘)  =  W  , 
W2  -W,  =  Wxs  . 

Then  equation  (3)  becomes 


it-  ptJpc  '-1  aisfi +#T[f;l0l 

-  —la  J 1  LAXaoJ 


(with  C  =  X*2  -  ,  according  to  Rummel  et  al.  (1978),  the  cross-track 

velocity)  or  as  a  formal  matrix  equation 


I=BY  +  A6£+«_ 


where 


rpc"iciT  T  r a Yioi  rp-lcV  T 

r  =  t  -  £  ,  B  =  “  ?  ,  Y  =  j  — 1  ,  and  A  =  ~  W  *  . 

L  J  LMaoJ  ~  J 

If  D  is  the  a  priori  variance  covariance  matrix  of  r  and  C  a  the  second-order 
moment  ("covariance  matrix",  reproducing  kernel)  of  with  6J  assumed  to 
be  an  element  of  an  infinite  dimensional  Hilbert  space,  then  the  best  linear  es¬ 
timates  for  Y  and  Jig  are,  seee.g.  Moritz  (1972): 

r&oi  [  rpa*ici  Y  rpe-lci 

H&JM  e„H  ejj  H  e,a  J  *  (6» 


M  •  Cg  WtJ 


9  -  as,  A’*D  =  [Peii£]  +£ 


Problems  in  using  this  type  of  least  squares  collocation  solution  come  from 
the  difficulties  in  deriving  the  matrix  of  equation  (8).  It  requires  the  connec¬ 
tion  of  the  mentioned  Hilbert  space  with  reproducing  kernel  (with  Cg  based  on 
a  chosen  degree  variance  model)  with  time  dependent  functionals  such  as  the  state 
vectors  of  the  two  satellites.  This  conceptual  problem  and  linked  to  it  the  problem 
of  the  practical  computation  of  Q  for  instance  by  numerical  or  analytical  inte¬ 
gration  will  not  be  investigated  in  this  report,  instead  an  alternative  approach  will 
be  analyzed  which  comes  also  closer  to  the  concept  of  "direct  gravity  mapping" 
as  formulated  by  Muller  &  Sjogren  (1968). 


3.  Observational  Model  for  Range  Rate  Changes  and  Orbit  Requirements 


Instead  of  working  with  the  range  rates  p  one  can  also  analyze  their  time 
derivative  p  ,  the  range  rate  changes.  Equation  (1)  yields; 

P(Si«Sa)  =  Xia  *  +  ^12  *  ®.ia  (9) 

For  SST  in  the  high-low  mode  Hajela  (1978,  p.  6)  showed  that  the  second  term  on 
the  right-hand  side  is  negligible,  a  fact  that  tremendously  facilitates  the  whole 
model.  It  is  to  be  shown  that  this  is  at  least  on  an  average  the  case  for  the  low- 
low  version,  too.  With 

2S is  '  £ia  =  P13  (  ~  &  )  ~  Pi2  C  . 

equation  (9)  becomes 

P(S\,Sa)  =  Xva  •  e,ia  +  P\a  Q3  (10) 

Introducing  the  reference  range  rate  change 

Pi 2  ~  ^£ia  *  ®ia  +  Pta  2 

as  derived  from  some  approximate  state  vectors  and  an  ellipsoidal  reference  field 
or,  in  addition,  a  set  of  known  potential  coefficients  and  assuming  e^s  and  pis 
to  be  known  (which  has  no  consequence  on  the  following  evaluations)  one  obtains 
as  difference; 

6p  =  5Xia  *  £ia  +Pa5£2-  <U) 

Evaluation  of  the  average  magnitude  of  the  second  term  on  the  right  hand  side  of 
equation  (11)  or  of  p“l  5xia  ,  since 

P»  <  Pi'a1  « 
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From  a  linearization  of  the  energy  conservation  law  follows  according  to  Wolff 
(1969): 

Mia  =  (Ma  -  Mi)2  -  T®  , 

where  X°  is  the  mean  velocity  of  the  satellites  and  T12  =  T  ( S2)  -  T  ( Sj)  the 
disturbing  potential  difference.  The  average  behavior  of  the  earth's  gravity 
potential  is  expressed  by  a  chosen  degree  variance  model,  with  o^(  T)  the 
disturbing  potential  degree  variances.  On  an  average  the  second  term  becomes 
therefore: 

ave<Pra  S*w>  "  ave(  T3)  =  - h-j  (var  T fo)  + 

Pis*  Pl3ii 

-  2  cov  (T(S1),T(Sa))  +  varT(S,)).  (12) 

Evaluation  of  the  average  magnitude  of  the  first  term  on  the  right  hand  side  of 
equation  (11),  i.e.  of 

6Xia  •  e  u  : 

E.g.  for  radial  separation  one  obtains 
...  a 


Mia*  *«  =  77(T(S  )  -T(Sx))  =  T,13  , 


ave(Mi2  *  «*)  *  (var  Tr  <si)  -  2cov(Tr(S1),  Tr(Sg))  +  var  Tr(Sa))  (13) 

The  average  magnitude  of  the  two  terms  on  the  right  hand  side  of  equation  (11)  for 
the  high  low  as  well  as  for  the  low  low  mode  is  contained  in  Table  1.  The  values 
have  been  determined  -  based  on  equations  (12)  and  (13)  -  from  equations  (A5)  and 
(A6),  described  in  the  appendix. 


Table  1.  Average  magnitude  of  the  two  terms  on  the  right  hand 
side  of  equation  (11) 


mode 

wave  lengths 

l 

first  term  of  (11) 
10” 5  msec”3 

second  term  of  (11) 
lO^msec”3 

high 

2  -  ® 

18.068 

0.277  •  10-6 

-low 

19  -  ® 

5.066 

0.000,*  10-6 

low 

2  -  ® 

1.498 

0.244  •  10”6 

-low 

19  -  « 

1.078 

0.016  •  10”6 

r 


The  numbers  of  Table  1  as  well  as  the  level  of  the  observational  noise, 
to  be  discussed  later,  clearly  indicate  that  the  second  term  can  be  neglected 
without  significant  loss  of  accuracy.  The  basic  model  will  therefore  be: 

P(SX,S8)  *  Xia  *  eia  (14) 

From  now  on,  range  rate  changes  are  considered  as  being  observed.  They  are 
obtained  from  range  rates  by  numerical  differentiation,  as  discussed  e.g.  in 
Rummel  et  al.  (1976)  and  Hajela  (1973).  The  mathematical  model  of  equation  (14) 
is  linked  to  the  observed  range  rate  changes  by 

=  P( ti;Xx<£).  Xa<£>,  Xi<£>.  Xs(£))  +  e,  (15) 

Assuming  for  a  moment  the  positions  Sx  and  S8  of  the  two  satellites  to  be  known, 
the  dependence  on  the  unknown  gravity  parameters  =  £  -  _5C  is  expressed  by 

f  i  —  Pi  (li  •  Si  ,S8)  +  6x  1 8  *  e^8  * 

and  since 

Hia  =  5Xa  -  6Xi  =  (gradiaT)(Sa)  -  (gradiiT)(S:) 

=  grad  T18  , 

the  observation  equation  becomes  in  vector  form: 

l  =  pc(t, ;  SltS8)  +  (gradg^(grad  T18))T  (S1 ,  S ^  6£  •  e ^  +  e  .  (16) 

In  analogy  to  the  last  chapter  the  uncertainty  of  the  satellite  positions  is  taken 
care  of  by 

i  =  Pc( t»;  St, S')  +  (grad,  p)T  (St)  ^Xj  +  (grad,  p)T  (S')dX3 

—1  *  * 

+  (gradgg  (grad  Ta  )r  (S' ,  S8)  65  •  e^  +  I  (17) 

The  terms  (grad^  p )T  (St )  and  (grad^  p)7(S3)  express  on  one  hand  the  coeffi¬ 
cients  for  an  orbital  improvement  from  the  observed  range  rate  changes,  on  the 
other  hand  they  define  the  sensitivity  of  the  observations  on  orbit  errors  as  ex¬ 
pressed  by  AX,  the  difference  between  the  actual  and  the  reference  orbit.  This 
type  of  coefficient  matrix  is  denoted  fl  ; 

O  =  (gradj^  P )  (Si)  =  (“  P  (Xia  “Pc  is)  +  [ 1  ei2  1  i,  j  =  1, 2, 3 

Expressing  fl  in  a  local  ( Z  >  £  ,  ‘H )  -  system  as  defined  in  Reed  (1973)  and  assuming 
radial  separation  of  the  two  satellites,  it  becomes 


>  ; 

t  ! 

::  I 
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n  -  -pc"  (Wu  -p  c  o  j) 


+  I  v,=^  i  <  -o' 


+  v,i 


where  e.g.  V12  stands  for  - 


is  now  split  into  the  contribution  of  the  normal  field,  fia  and  into  the  average 
contribution  of  the  residual  field,  Qr  ,  i.e.  : 

n  =  n «  +  Hr  and 


=  ave  (  -  p 


r 

1  12 


/ 

)  +  ave  T £* 


Magnitude  of 

The  T-terms  for  the  determination  of  fiB  are  (Heiskanen  and  Moritz,  1967, 
p.  231): 


rC  =  Z*(2l+  1)  J3£(!)a*Pai(sin0)  , 

r  £=i  v r/ 

r«  .  gm  ra  /ay3*  dp3* 

r  ^  ,  =!  Ja*  \r/  d®  ’ 


rn  =  o , 


with  GM 

Ja£ 

a 


GM  . . .  geocentric  gravitational  constant 

Ja£  ...  spherical  harmonic  coefficients 

a  ...  semimajor  axis  of  the  earth 

P2£  (cos  ip) . . .  Legendre  polynomial. 


Inserting  the  geodetic  constants  of  the  Geodetic  Reference  System,  1967  and 
assuming  a  radial  separation  of  50  km  at  200  km  altitude  one  obtains  at  most : 


j  -0.272  •  10 "6  sec-2  radial  (p  =  0°) 

£2a  =  j  0.000  84  •  10"5  sec-2  latitudinal  (<jfl  =  45°) 

(  0  longitudinal 

Foranalongor  cross  track  separation  all  three  components  are  less  than 
1  E.U.  =  10~9  sec-2  which  is  obvious  because  of  the  distinct  dominance  of  the 
radial  term  in  the  reference  field. 


Magnitude  of  fir 

r  rz 

The  average  values  for  T  ^  and  are  again  taken  from  an  adopted 

degree  variance  model.  Equations  (A6)  and  (A7)  of  the  appendix  yield 

C  r  f  j* 

=  ave  (-p_1  Tu)  +  ave  T12  =  0.63  E.U. 

The  two  other  components  will  on  the  average  be  at  most  of  the  same  magnitude. 

The  entire  contribution,  6 p,  of  an  orbit  error  &Xx  or  A_xx  if  expressed 
in  the  (C  »  £  .  V  )~  system  is; 

6o  =  (f2n  +  nr)7  Ax  . 

The  orbit  corrections  AXi  and  AXa  depend  on  the  gravity  parameters  60  (short 
wavelength  part)  and  on  the  initial  state  vector  corrections  AXl0  and  AX*  (long 
wavelength  part)  and  so  does  Ax  .  Of  concern  for  the  gravity  parameter  estimation 
are  unmodelled  short  wavelength  effects.  They  are  caused  especially  by  non- 
gravitational  disturbances,  and  to  a  less  extent  (because  of  the  long-wavelength 
characteristic)  by  tracking  error  and  tracking  station  uncertainties.  The  latter 
when  multiplied  with  short -wavelength  coefficients,  i.e.  with  Or  ,  may  cause 
problems,  too.  In  figure  1  the  magnitude  of  the  range  rate  change  error  due  to 
an  orbit  error  for  the  terms  (radial),  (along  or  cross  track),  and 
n£(=  nr£=  a  Jl)  are  shown.  It  clearly  indicates  that  especially  the  disturbances 
in  the  radial  component  of  the  orbit  are  critical.  There  are  also  given  three  levels 
of  measurement  precision  (0.001,  0.01,  and  0. 1  mgal).  Their  relation  to  what  is 
actually  feasible  will  be  given  later. 
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0.01  0.1  1.0  10.0  AX  m 

Figure  1.  Range  rate  change  error  due  to  an  orbit  error  for  the  terms 

£2>  (radial),  flaS  (along  or  cross  track),  and  o£(s=!n,£*!  JW 


In  the  following  -  after  having  some  insight  into  the  impact  of  errors  on 
the  orbit  -  the  orbital  correction  will  be  disregarded.  From  equation  (17)  the 
estimation  model  becomes  therefore: 

1  =  2e(ti;SlfSa)  +  (grad^g  (grad  Tm))t  (St ,  S3 )  +  e  ( 


or  with 

£  =  i  -Oc(tt  jSi.Sjj),  and  A  =  (gradgg  (grad  T^))1  (S;,  S2): 

r  =  A  6ff  +  e  (19 

With  D  now  the  a  priori  variance  covariance  matrix  of  the  range  rate  changes 
r  and  Cg  as  before  the  best  linear  estimate  of 


min  (  |1  r  -  A  +■  |'|  Sg'l^-i  ) . 


becomes 


M  *  (AT  D"1  A  +  C"1)"1  At  D"1  r 
=  C&Ar  (A  C.At  +  D)~l  r 


=  L  r 


where  L  =  Cg  Ar  (A  CgAJ  +  D)~l  is  the  estimator. 


4.  Computational  Procedure  and  Error  Measure 


The  goal  of  this  study  is  to  perform  an  error  analysis  of  SST  in  the  low-low 
mode.  The  adequate  error  measure  is  the  a  posteriori  variance-covariance  matrix 
for  the  gravity  parameters  to  be  estimated  from  a  realistical  arrangement  of 
sample  points.  The  a  posteriori  variance-covariance  matrix  provides  on  one  hand 
a  reliable  measure  of  the  variance  of  the  quantities  to  be  estimated  and  on  the 
other  hand  the  covariance  structure,  too,  which  shows  how  independent  the  esti¬ 
mates  actually  are.  The  commonly  used  form  of  the  a  posteriori  variance- 
covariance  matrix,  derived  from  equation  (21),  is: 

Eg  ~  —8  ~  £<9  — T  A  T  +  D )  1  A  Cg  .  (22) 

A  A 

Eg  is  the  second-order  moment  E  f  [  <3/3  -  6£]  [6g  -  6£]  } .  With  model  equation 
(19)  and  the  estimate  (21)  one  finds 

Eflfig- =  E[[L  (A  6£  +  e)  -  8gHL  (A  6g+e)  ~  *glT} 

■  B  E{5£  5gT}  Bt  +  LDLt  , 

where  B  s  (1,  A  -  J)  ,  ,  identity  matrix,  and  E  {  e  }  =  0.  The  expectation  opera¬ 

tor  has  to  be  replaced  by  the  averaging  operator,  M,  over  the  sphere,  see  e.g. 
Heiskanenand  Moritz  (1967,  p.258)  to  obtain 
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(23) 


( 

! 


Eg  =  B  Cg  BT  +  LD  Lt 

This  is  an  alternative  formulation  of  equation  (22).  It  presents  with  the  first 
term,  B  B  \  the  configuration  inadequacy;  B  =  LA  -_I  would  be  0  in  the 
optimal  case.  The  second  term  shows  the  propagation  of  the  measurement  error. 
Thus  the  contribution  of  the  first  term  could  be  minimized  by  optimizing  the  con¬ 
figuration,  that  of  the  second  term  mainly  by  improving  the  observational  pre¬ 
cision. 


In  equation  (22)  A  Cg  Ar  is  nothing  but  C  ^ ,  the  a  priori  variance  covar¬ 
iance  matrix  of  the  observations,  and  Cg  At  the  cross -covariance  Cg-p  between 
the  parameters  to  be  estimated  and  the  range  rate  changes,  or: 

Efl  -  Cp-CpfiiCp  +  D)_1C^  (24) 

Since  from  equation  (14) 


P  =  <Xa  -  X^-e^  =  (7T2  -  vTl)-eu 

where  T  will  be  referred  to  a  known  reference  field  up  to  degree  ^nax  =  12,  a  single 
matrix  element,  [ijj,  of  Cp  becomes 

CjM  =  ei,iC[(VTa  -  VTx)t;  ( 7T3  -  7Tt)j  }  £u,j 

=  eja.t  [CfVT,I,;7Xl>,}-Cf7T.,,;7T1#J}  -  CfvTlf ,  ;7Ts,j} 

+  C{VTv;7T1(J)]el2)J  .  (25) 

The  gradient  is  given  in  an  earth-fixed  (  r ,  <p  ,  A)  -  system  with 

3T  1  3T  1  3T  1T 

dr  *  r  d<fi  ’  r  cos  «)  5A  ! 

Thus,  before  inserted  into  (25)  the  unit  vectors  eu  provided  in  an  earth-fixed 
(x,  y,  z)  -  system  have  to  be  transformed  to  the  (  r,  <fi  ,  A)  -  system  (see  Hajela, 
1978,  p.  8)  by 


sin  ip~ 
cos  <0 


^  (x,  y  ,  z) 

12  _  _13 


£52 


and 


cos  <0  cos  A 
-  sin  0  cos  A 

-  sin  A 


cos  <p  s  in  A 
s  in  O  s  in  A 
cos  A 


Similarily  an  element  of  becomes: 


[CC«0,7TafJ]  -  Cf6|3,7TltJ}lr  eUfJ  . 


(26) 
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The  subroutine  COVAX  described  in  Tscheming  (1976)  provides  all  needed 

/>T  1  vr  1  VI*  \ 

covariance  elements.  Instead  of  the  gradient  elements  —  •-  ~  coa  ~  7^  ; 

1  >T  1  5T  1  5T 

e.g.  In  units  of  mgal  [  =  1<T  msec  ].  it  gives  “  •  ”77  "'yr  Cos«  TT  / 

in  units  of(E.U.t  arcsec,  arcseci,  with  y...  normal  gravity.  The  conversion  is 
carried  out  by  means  of  the  diagonal  matrix 


-  r  10* 
0 
0 


-<y  o“) 


0 

0 

-o/o* 


The  underlying  degree  variance  model  used  in  all  computations  is  isee  Tscherning 
and  Rapp,  1974) : 


_  425.28  (£-1 


[  mgal2 


( i  -  2)(  £  +  24) 

for  the  gravity  anomaly  degree  variance,  and 
r.  =  6. 369  779  8  •  10®  m 


for  the  radius  of  the  Bjerhammar  sphere.  The  computation  of  one  single  element 
of  Cp  requires  6x4  different  covariance  elements  and  is  therefore  very  time 
consuming. 

As  an  alternative  we  used  therefore  the  subroutine  COVAPP  of  Siinkel  (1979). 
This  subroutine  computes  the  covariance  elements  by  a  spline  interpolation  from 
elements  stored  once  for  all  in  a  table.  The  access  is  organized  almost  identical 
to  COVAX.  For  a  sufficiently  dense  interpolation  table  both  subroutines  yield 
almost  identical  results.  COVAPP  is  much  faster  than  COVAX  but  requires,  of 
course,  additional  core  storage. 

As  described  in  the  introduction  the  gravity  parameters  investigated  in  the 
error  analysis  will  be 

-  geoid  heights  and  geoid  height  differences,  and 

-  l°x  1°  mean  gravity  anomalies  (which  contain,  roughly,  information 
down  to  a  half-wavelength  of  100  km). 

The  estimates  for  mean  gravity  anomalies  are  obtained  by  averaging  over  81 
point  estimates  in  each  l°x  1°  block. 


It  remains  to  define  the  a  priori  variance  covariance  matrix  D  of  the 
observations.  In  a  real  world  experiment  one  could  depart  from  the  original  range 
rate  observations  and  derive  from  them  range  rate  changes  by  numerical  differen¬ 
tiation,  as  described  in  Rapp  and  Hajela  (1979).  Since  we  are  going  to  work  directly 
with  range  rate  changes  we  have  to  define  precision  levels  for  them  which  cor¬ 
respond  in  a  well  defined  manner  to  range  rate  precisions.  We  define  that  certain 
range  rate  and  range  rate  precisions  correspond  to  each  other  if  the  same  piece 


of  information,  e.g.  mean  gravity  anomalies  with  approximately  the  same 
standard  deviation  and  the  same  resolution  can  be  derived  from  them.  In  order 
to  deduce  the  equivalent  precision  levels  the  degree  variance  model  of  one 
observed  quantity,  e.  g.  range  rate  changes  equation  (A6),  is  compared  with 
the  corresponding  degree  variance  model  of  the  observation  noise.  The  resolu¬ 
tion  of  the  measurement  process  is  then  defined  by  that  frequency  at  which  the 
signal,  represented  by  the  degree  variance  model,  is  equal  to  the  noise  in 
amplitude.  For  higher  frequencies  the  noise  would  dominate  the  signal.  In 
order  to  have  the  same  resolution  for  the  observation  type  to  be  compared  with, 
e.g.  range  rates,  the  signal-to-noise  ratio  has  also  to  be  equal  to  one  at 
this  "resolution  frequency".  Under  this  assumption  one  can  deduce  from  the 
degree  variance  model  of  the  signal  and  noise  of  the  second  process  -  equations 
(AS)  and  (A9)  -  fee  corresponding  precision  of  the  second  process.  Based  on 
the  expressions  for  the  degree  variance  models  of  signal  and  noise  for  range 
rate  and  range  rate  change  -  derived  in  fee  Appendix  -  figure  2  is  computed. 

The  graph  relates  range  rate  and  range  rate  change  ( *  acceleration  in  the  line 
of  sight)  precisions,  covers  a  range  from  an  altitude  of  150  km  to  300  km  of 
fee  satellite  system,  and  is  valid  for  low-low  and  high-low  configurations. 


icr9  XT8  icr7  10- 6  10'5  10H  ms'1 

Figure  2.  Range  rate  change  (acceleration  in  the  line  of  sight)  precision 
versus  range  rate  precision. 


Figure  2  enables  us  to  assign  to  each  range  rate  precision  introduced  into 
the  a  priori  variance  covariance  matrix  D  the  equivalent  precision  for  observed 
range  rates.  The  observations  are  assumed  to  be  uncorrelated  and  of  equal 
variance  a*  ,  or  D  =  ad _I .  Representative  range  rate  precisions  which  are 
envisaged  for  coming  experiments  and  the  equivalent  precision  for  range  rate 
changes  are  e.g.  : 

10“*  ms'1  ■  1.5  •  10-6  ms“a 
10“*  ms-1  =  2.2  •  10 “®  ms"2 
10"*  ms-1  =  3.1  •  10"’’  ms'3 
10'7  ms-1  =  4.0  •  10 ms"2 

These  numbers  allow  now  to  valuate  the  magnitude  of  the  second  term  of  equation 
(11)  (see  table  1),  and  the  impact  of  systematic  and  random  errors  in  the  orbit 
as  displayed  in  figure  1. 

For  numerical  simplicity  Cjj  in  equation  (24)  is  decomposed  into  its 
eigenvalues: 

Cp-  =  U  A  Ur  (2',) 

(U...  matrix  if  eigenvectors  of  C-  with  UV*  =UU  =_I ,  and  A  ...  diagonal 
matrix  containing  the  eigenvalues  A,).  As  discussed  in  Rummel  et  al.  (1979)  this 
allows  on  one  hand  the  repetition  of  the  error  analysis  with  various  observation 
noise  levels  <70a  with  solving  the  system  of  linear  equations  (27)  only  once  and 
gives  on  the  other  hand  a  perfect  insight  into  its  stability  behavior. 

Equation  (24)  becomes  then: 

Eg  -  C0  -  Cg-p  u  (A  +  ac3!)-1  u  Cgja 

or  with  Cg^  y  =  V: 

Eg  =  -  V  (A  +  (703  Ij"1  Vr.  (28) 

where  (A  +  a,2^)"1  is  a  diagonal  matrix  with  elements  l/(At  o02).  Fora 
separation  of  the  contributions  of  BCgB:  and  LD  L1  of  Eg  in  equation  (23) 
the  decomposition  is  also  performed  for  LD  LT  : 

L  D  Lr  =  £ SH  ( — p  +  <?o3I  <C6-  +  <7oa±)~l  Cfo 

-  v  (A  +  a03_i)_1  a03I(A  +  Oo3!)"1  vT  .  (29) 

Finally,  it  is  necessary  to  generate  sample  points  for  the  two  satellites. 

A  circular  orbit  has  been  generated  with  the  following  characteristics: 


mission  period  50  days 

inclination  68° 

sampling  interval  20  sec 

altitude  of  S.  200  km  (and  150  km,  250  km,  300  km) 

period  5301  sec 


For  stability  reasons  spherical  distances  of  smaller  than  0.3°  between  two 
adjacent  sample  points  were  not  permitted.  The  area  for  which  sample  points 
are  used  in  the  analysis  is  (X  |  0°  =»  X  s  5°  ]  and  f«J  |  1°  *  o  <  9°  ]  .  The 
resulting  ground  track  coverage  is  given  in  figure  3. 

Figure  3  also  contains  the  six  l°x  1°  blocks  for  which  the  mean  grav  ity 
anomaly  error  analysis  will  be  performed.  In  order  to  pronounce  the  special 
features  of  a  radial,  along  track,  or  cross  track  arrangement  of  the  two  satel¬ 
lites  the  second  satellite  is  kept  exactly  in  one  of  these  modes  throughout  the 
entire  mission  period  for  each  run  to  be  analyzed. 


5.  Results 


All  a  posteriori  standard  deviations  for  mean  gravity  anomalies,  geoid 
heights,  and  geoid  height  differences  will  refer  to  a  reference  field  of  degree 
and  order  12,  assumed  to  be  known. 


l°x  1°  Mean  Gravity  Anomalies 

Figure  4  shows  the  standard  deviation  of  an  estimated  l°x  1°  mean  gravity 
anomaly  (block  no.  2  in  figure  2)  computed  from  the  74  sample  points  given  in 
figure  2  as  a  function  of  the  assumed  range  rate  change  observational  precision 
<7(P) .  The  corresponding  range  rate  standard  deviation  is  also  indicated.  The 
graphs  are  given  for  a  cross  track,  along  track,  and  radial  separation  of  250  km 
of  the  two  satellites  at  analtitude  of  200  km. 


The  attainable  accuracy  of  ±7.5  mgal  for  a  radial  separation  is  much  better 
than  that  derived  from  an  along  track  separation  which  is  again  better  than  that  of 
a  cross  track  separation.  At  a  first  glance  this  result  surprises.  The  structure 
of  the  gravity  field  in  radial  direction  should  be  not  much  different  from  that  in 
horizontal  direction.  For  instance  the  degree  variance  <r/(  6r)  of  the  radial 
component  of  the  gradient  of  the  disturbing  potential  relates  to  that  of  the  horizon¬ 
tal  components  ct^3(6^)  or  cra(6,)  (longitudinal  and  transversal),  see  equations 
(A3)  and  (A 4)  of  the  appendix: 


0^(5r)  =  2 


l  +  l 
l 


r°t(h) 

S3<6-) 


for  l  >  12. 


The  contribution  of  the  pure  error  propagation  part  -LDL1  of  equation  (23)  - 
also  contained  in  figure  4  is  the  same  for  all  three  types  of  separation  and  almost 
as  low  as  the  theoretical  lower  bound  discussed  in  Rummel  (1978).  Thus  the  worse 
results  for  along  and  cross  track  have  to  be  caused  by  a  configuration  problem. 

The  correlation  functions  for  acceleration  differences  at  an  altitude  of  200  km  - 
figures  5a  and  5b  -  and  for  the  cross -covariance  between  a  surface  gravity  anomaly 
and  acceleration  differences  -  figures  6a  and  6b  -,  both  for  separations  of  10  km 
and  of  250  km  provide  an  answer:  The  area  covered  by  observations  -  see  figure 
2  -  has  only  an  extension  of  8°x  5°.  But  whereas  the  correlation  functions  for 
radial  separation  have  a  zero  at  approximately  5°  and  little  power  at  spherical 
distances  greater  than  5°  this  is  not  at  all  the  case  for  along  track  separation  and 
even  worse  for  cross  track  separation.  Thus,  the  spatial  arrangement  of  the  two 
satellites  does  not  principally  influence  the  results  of  the  gravity  parameter  esti¬ 
mation.  For  an  along  track  or  cross  track  arrangement  only  a  larger  area  around 
the  estimation  points  has  to  be  covered  by  observations  to  yield  results  of  the 
same  quality  as  obtained  for  radial  separation. 


Figure  5a.  Covariance  function  of  acceleration  differences  for  a 
satellite  altitude  of  200  km  and  a  separation  of  250  km 
(Pi  with  0  =  0°  and  A=0°;  Pa  with  =  and  A=0°) 


Figure  5b.  Covariance  1 
satellite  alti 
(PL  with  <p  = 


'unction  of  acceleration  differences  for  a 
tude  of  200  km  and  a  separation  of  10  km 
0°  and  X  =  0°  ;  P3  w  ith  =  0  and  X  =  0° ) . 


Figure  6a.  Cross-covariance  function  between  a  surface  gravity  anomaly  and 
acceleration  differences  at  an  altitude  of  200  km  and  a  seperation 
of  250  km  (P  with  «J  *  08  and  X  =  0°;  Pa  with  o  =  ij>  and  X  =  0°) 
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Figure  6b.  Cross-covariance  function  between  a  surface  gravity  anomaly  and 

acceleration  differences  at  an  altitude  of  200  km  and  a  separation  of 
10  km  (P:  with  <D  =  0°  and  A  =0°;  P3  with  <fi  =  0  and  A  =0°). 


Next,  the  impact  of  the  separation  of  the  two  satellites  on  the  standard 
deviation  of  the  estimated  mean  gravity  anomaly  has  been  investigated.  Figure 
7  contains  the  graphs  for  (T(^g)  for  a  radial  separation  of  10  km,  50  km,  and 
250  km.  In  all  cases  the  lower  satellite  is  at  an  altitude  of  200  km.  The  contri¬ 
bution  of  LD  l/  ,  equation  (23),  is  given,  too.  As  to  be  expected,  the  accuracy 
of  the  estimated  gravity  anomalies  improves  with  increasing  separation,  if 
CT(p),  or  cr((5)  respectively,  is  kept  constant.  This  result  is  discussed  already 
in  Rummel  (1978)  based  on  purely  theoretical  considerations. 
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But  not  only  the  standard  deviation  of  the  estimates  should  be  analyzed. 
Small  correlations  are  a  measure  of  a  high  mathematical  independence  of  the 
estimated  parameters.  For  block  no.  2  (figure  2)  and  for  an  optimum  of  obser¬ 
vational  precision  attainable  in  a  gravitational  satellite  mission  -  ±0.03  mgal 
for  a(p)  or  ±  10"°  ms-1  for  <7(  p  ),  respectively  -  the  standard  deviation  of  the 
a  posteriori  estimates  and  the  correlations  are  presented  in  tables  2a  through  2f. 
They  contain  on  the  main  diagonal  the  standard  deviation  of  all  six  considered 
1x1°  blocks  and  on  the  off-diagonals  the  correlations  among  them.  The  results 
are  given  for  radial,  along  track,  and  cross  track  separations  of  10  km  and  250 
km  at  an  altitude  of  200  km.  For  radial  separation  almost  perfect  independence 
is  achieved.  The  results  for  along  track  and  cross  track  separation  show  a 
somewhat  higher  correlation  in  the  direction  of  the  separation.  Again,  this 
higher  correlation  would  decrease  with  data  covering  a  larger  area. 


Table  2a.  Standard  deviation  and  correlations  for  the  six 
l°x  1°  mean  gravity  anomalies. 


1 

2 

3 

4 

5 

6 

8.4 

-Q23 

0.11 

i 

.o 

-0.10 

0.06 

1 

8.2 

-0.23 

1 

O 

o 

1 

O 

-0.11 

2 

8.3 

0.06 

0.11 

0.14 

3 

8.4 

OO 

CM 

O 

_ l _ 

0.13 

4 

radial 

<r(  p)  =  t  0 . 03  mgal 
( tz  a(  p)  =  -  10"6  ms-1 ) 

8-1 

026 

5 

rx  =  200  km 
sep.  10  km 

8.2 

6 
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Table  2b.  Standard  deviation  and  correlations  for  the  six 
l°x  1°  mean  gravity  anomalies. 


1  2  3  4  5  6 


12.7  -0.02-0.03  0.35  0.02  0  1 


13.3  0.10  023  0.41  0.02  2 
13.5  0.08  0.27  0.40  3 


along  track 
cr(p)  =  ±  0.03  mgal 
(  =  a(P)  -*  10 “6  ms-1) 
rt  =  200  km 
sep.  10  km 


12.8  0.07  0.05 


13.1  0.11  5 

13.4  6 


Table  2c.  Standard  deviation  and  correlations  for  the  six 
l°x  1°  mean  gravity  anomalies. 


1 

2 

3 

15.1 

Q55 

0.52 

cross  track 
c r(p)  =  1  0.03  mgal 
(icr(f&)  =i  10"a  mi 
rt  =  200  km 
sep.  10  km 


5  6 


0.34  1 


0.53  0.26  0.19  0.25  2 


14.8  0.34  0.25  0.27  3 


15.3  0.52  0.52 


14.6  0.531  5 


14.8  6 


1 

2 

3 

a 

5 

6 

6.9 

0.09 
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Table  2e.  Standard  deviations  and  correlations  for  the  six 
l°x  1°  mean  gravity  anomalies. 


1 

2 

3 

D 

5 

6 

9.3 

0.12 

0.34 

0.04 

015 

0.26 

1 

9.9 

0.02 

0.27 

0.07 

-0.01 

2 

10.4 

0.19 

0.16 

009 

3 

10.4 

0.11 

0.13 

4 

along  track 

j(p)  *  i  0.03  mgal 

[  £  cr<  6)  *  *  io"8  ms'1) 

10,8 

0.02 

5 

rx  ■  200  km 
sep.  250  km 

11.3 

6 

-26- 


Table  2f.  Standard  deviations  and  correlations  for  the  sue 
l°x  1°  mean  gravity  anomalies. 
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Another  important  factor  influencing  the  quality  of  the  estimated  gravity 
parameters  is  the  altitude  at  which  the  experiment  is  performed.  Figure  8  gives 
the  results  analogous  to  figures  4  and  7  but  now  with  a  constant  radial  separation 
of  50  km  and  for  altitudes  of  150  km,  200  km,  250  km  and  300  km.  The  attainable 
±5  mgal  for  a  l°x  1°  mean  gravity  anomaly  from  only  74  observations  of  an  al¬ 
titude  of  150  km  would  make  it  definitely  worthwhile  to  fly  the  experiment  at  this 
low  altitude  if  the  related  technical  problems  could  be  controlled.  The  increase 
of  the  accuracy  with  decreasing  altitude  is  very  distinct.  Table  3  contains  the 
results  for  all  six  blocks  analogous  to  tables  2  for  an  altitude  150  km. 
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Figure  8.  Variation  of  the  estimated  a  posteriori  standard  deviation 
of  a  l°x  1°  mean  gravity  anomaly  with  changing  experiment 
^  .  altitude  (radial  separation  of  50  km)  and  the  contribution  of 

IQ]  L  D  Lr  (eqn.  (23)). 


0.1  a  (§)  1.0  10.0  *10'6ms'1 


Table  3.  Standard  deviations  and  correlations  for  the  six 

l°x  1°  mean  anomalies  for  an  experiment  altitude 
of  150  km  and  a  radial  separation  of  50  km. 


Finally,  with  a  less  optimal  arrangement  of  the  sample  points  with  respect  to 
the  i°x  1°  blocks  -  shown  in  figure  9  -  the  whole  analysis  has  been  repeated. 
Representative  for  these  results  are  the  numbers  for  the  six  blocks  given  in 
table  4  derived  for  a  radial  separation  of  250  km  at  an  altitude  of  200  km.  They 
show,  when  compared  with  table  2d,  a  significant  decrease  in  the  accuracies  of 
approximately  20%. 
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Table  4.  Standard  deviations  and  correlations  for  the  six  ^ 

l°x  1°  mean  anomalies  with  a  suboptimal  -j 

arrangement  of  the  sample  points.  I  •' 
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Geold  Heights  and  Geoid  Height  Differences 

In  figure  10  the  ground  tracks  are  shown  together  with  the  locations  for 
which  the  a  posteriori  variances  of  geoid  heights  and  geoid  height  differences 
have  been  estimated  (x)  .  Also  contained  are  so  called  suboptimal  locations 


mm 


0.0°  1.0°  2.0°  3.0°  4.0°  5.0°  X 

Figure  10.  Groundtracks  with  the  locations  of  the  points  for  which 
geoid  heights  are  estimated  (”x”  optimal  locations; 

"A”  arbitrary  arrangement). 


i 


The  estimation  of  geoid  height  differences  is  discussed  in  Christodoulidis  (1976) 
and  Jekeli  (1979).  The  a  posteriori  variance  of  a  geoid  height  difference  6n 
between  two  points  Pt  and  P3  is  derived  from  the  variances  and  covariances 
of  the  point  estimates  by: 

<73  (6N(P3  -  Pr)}  =  <7s[N(Pa)}  -  2  cov{N(Pa);N(P1)}+a2fN(P1)}  (30) 

Equation  (30)  indicates  that  mainly  high  correlations  are  responsible  for  low 
variances  or  standard  deviations  of  the  geoid  height  differences.  To  me  it  seems 
therefore  no  principal  advantage  to  estimate  differences  of  a  quantity  instead  of 
the  quantity  itself,  except  for  the  aspect  of  eliminating  systematic  errors  with 
long  wavelength  structure.  Quite  in  contrast,  equation  (30)  shows  that  for  inde¬ 
pendent  estimates  the  standard  deviation  of  a  difference  is  /l  times  that  of  the 
point  estimate.  The  characteristic  of  estimating  differences  is  to  be  distinguished 
from  measuring  differences,  as  e.g.  in  gravimetry,  where  the  advantage  may  come 
from  the  measurement  principle. 

The  results  will  be  given  for  point  no.  3  and  for  geoid  height  differences  between 
points  no.  3  and  4  (151  km)  and  points  no.  3  and  5  (302  km).  Figure  11  shows  the 
a  posteriori  standard  deviation  of  the  estimated  geoid  height  N(P3)  and  geoid 
height  differences  6N(P4-P3)  and  6N(P3-P3)  as  a  function  of  the  observational 
precision  a  ( (3).  The  experiment  altitude  is  200  km  and  the  separation  250  km  in 
radial,  along  track  and  cross  track  direction. 

Again,  observational  data  covering  a  largerareawould  produce  estimates  of 
similar  quality  for  along  and  cross  track  separation  as  for  radial  separation.  The 
optimal  number  of  approximately  ±0.70  m  is  achieved  for  6N(  3-4)  in  radial 
separation.  As  to  be  expected  from  equation  (30)  the  results  for  geoid  height 
differences  depend  very  much  on  the  correlation  between  the  point  estimates.  This 
produces  usually  a  pattern  of  increasing  a  posteriori  standard  deviations  for  geoid 
height  differences  with  increasing  distance  of  the  two  points,  and  point  estimates 
slightly  worse  than  the  geoid  height  differences  for  small  distances.  Especially 
good  estimates  for  geoid  height  differences  are  obtained  when  the  two  points  are 
located  along  the  direction  of  separation  of  the  two  satellites,  either  along  or 
cross  track. 

The  dependence  on  the  separation  distance  of  the  two  satellites  is  displayed 
in  figure  12.  In  all  cases  a  radial  separation  of  an  altitude  of  200  km  is  assumed. 

As  for  mean  gravity  anomaly  estimates  the  results  improve  with  increasing  inter- 
satellite  distance. 
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Figure  12.  Variation  of  the  estimated  a  posteriori  standard  deviation 

of  the  geoid  height  N(P3)  and  the  geoid  height  differences 
/A\  ®N(P  -P3)and  6N(PS -  P3)  with  changing  separation  of  the 

O  IIN/^  UlOINy  two  satellites  (altitude:  200  km;  radial  separation). 
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Rather  heterogeneous  looks  the  picture  for  the  analysis  of  the  dependence 
of  the  estimates  on  the  experiment  altitude,  shown  in  figure  13.  The  standard 
deviation  of  the  point  estimates  increase  with  decreasing  experiment  altitude. 

The  reason  for  this  paradox  result  is  again  the  comparably  small  area  covered 
by  the  sampling  points.  For,  the  cross -covariance  function  between  surface 
geoid  heights  and  the  radial  component  of  the  disturbing  potential  shows  a  first 
zero  at  about  7°  with  a  comparably  high  correlation  for  large  spherical  distances. 
For  geoid  height  differences  the  a  posteriori  standard  deviation  would  agree  with 
the  anticipated  behavior  -  decreasing  standard  deviation  with  decreasing  experi¬ 
ment  altitude  -  but  only  for  very  low  range  rate  change  or  range  rate  noise  level. 

Finally,  the  results  for  points  located  along  a  ground  track  (  "  x  "  in  figure 
10)  are  compared  with  those  for  points  located  arbitrarily,  somewhere  between 
the  ground  tracks  (  "  A  "  in  figure  10).  Table  5  contains  the  a  posteriori  standard 
deviation  for  three  point  geoid  height  estimates  and  two  standard  deviations  for 
geoid  height  differences,  all  for  a  radial,  along  and  cross  track  separation  of 
250  km  at  200  km  altitude,  and  for  a  radial  separation  of  50  km  at  150  km  and  200 
km.  The  left  of  each  pair  of  columns  shows  the  numbers  for  the  optimal  point 
locations,  the  right  columns  those  for  the  arbitrary  locations.  A  significant 
dependence  on  the  point  location  can  only  be  seen  for  radial  separations  whereas 
for  along  and  cross  track  separations  the  point  lo nation  plays  no  distinct  role. 

Table  5.  Estimated  standard  deviation  for  geoid  heights  and  geoid  height 
differences  assuming  a  range  rate  change  precision  of  -0.03 
mgal  for  an  optimal  and  an  arbitrary  location  of  the  estimation 
points,  see  figure  10  (altitude/separation). 


Figure  13.  Variation  of  the  estimated  a  posteriori  standard 
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Information  Content  Per  Frequenc; 


Of  special  importance  for  the  user  is  that  the  estimated  gravity  parameters 
really  contain  information  down  to  a  half  wavelength  of  100  km,  or  up  to  l  =  180 
in  terms  of  spherical  harmonics.  We  computed  therefore  from  the  given  data 
the  variance  and  standard  deviation  of  the  geoid  height  of  point  no.  3  for  the  fre¬ 
quency  ranges  [l  1 13  £  £  <  “  },  { 2  |  121  £  l  <  00 },  and  { l  j  181  ^  l  <  * } 
separately  and  compared  the  information  contained  in  the  ranges  13-120,  121-180, 
181-“ .  Table  6a  gives  the  values  for  the  a  priori  standard  deviation,  a  posteriori 
standard  deviation  and  the  ratio  of  a  posteriori  to  a  priori  variance  for  the  three 
spectral  ranges  for  a  radial  separation  of  250  km  and  an  altitude  of  200  km.  Table 
6b  contains  the  corresponding  information  for  a  radial  separation  of  10  km. 

The  ratios  of  the  a  posteriori  to  a  priori  variance  indicate  that 

-  considerable  information  is  contained  in  the  frequency  range  from  121  to  180 

(1:15), 

-  almost  no  information  has  been  extracted  from  frequencies  higher  than  L  =  180 
(1:1.2  and  1:1.7  respectively) 

-  the  improvement  of  the  larger  separation  of  250  km  against  that  of  only  10  km 
comes  only  from  the  low  frequencies  (1:52  versus  1:22). 


Table  6a.  A  priori  and  a  posteriori  standard  deviation  and  ratio  of 

the  a  priori  to  the  a  posteriori  variances  for  three  spectral 
ranges  (altitude:  200  km;  radial  separation  250  km). 
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Table  6b.  A  priori  and  a  posterori  standard  deviation  and  ratio  of 

the  a  posteriori  to  the  a  priori  variances  for  three  spectral 
ranges  (altitude:  200  km;  radial  separation  10  km). 


spectre 

13-120 
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2S 
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1. 17 
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0. 42 
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1:22 
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6.  Summary  and  Conclusions 


The  goal  of  this  study  was  an  error  analysis  for  lc  x  1°  mean  gravity 
anomalies,  point  geoid  heights,  and  geoid  height  differences  derived  from  ob¬ 
servations  in  a  'low-low"  satellite-to-satellite  tracking  experiment.  The  chosen 
estimation  method  was  least  squares  collocation.  In  order  to  facilitate  the  mathe¬ 
matical  model  not  the  range  rates  themselves  but  their  time  derivatives,  i.e. 
range  rate  changes  or  line  of  sight  accelerations,  were  assumed  to  be  given. 

The  orbit  requirements  were  treated  separately  in  an  empirical  sensitivity 
study.  If  we  assume  a  range  rate  change  precision  of  =0.03  mgal  ( *  10  °  m  s'1 ) 
and  want  to  keep  systematic  errcrsbelow  1/10  of  this  level,  unmodelled  high  fre¬ 
quency  disturbances, in  the  orbit,  e.  g.  due  to  drag  effects,  would  have  to  be  kept 
below  1  cm  in  radial  direction,  and  1.  5  m  in  along  or  cross  track  direction. 

Thus,  a  drag-free  capability  or  a  system  of  micro-accelerometers  seems  to  be 
necessary.  Unmodelled  long  wavelength  errors  modulated  by  short  wavelength 
model  coefficients  have  to  be  kept  smaller  than  around  10  m  independent  of  the 
direction  of  the  separation,  and  should  therefore  be  no  limiting  factor. 

The  a  posteriori  standard  deviations  of  the  l°x  1°  mean  gravity  anomalies 
were  analyzed  as  a  function  of  the  experiment  altitude,  the  separation,  the  spatial 
arrangement  of  the  two  satellites,  and  the  location  of  blocks  with  respect  to  the 
sample  points.  For  a  separation  of  250  km  of  200  km  altitude  and  a  range  rate 
precision  of  ±0.03  mgal  a  standard  deviation  of  ±6.5  mgal  was  derived,  with 
almost  no  correlation  between  the  individual  blocks. 

The  same  criteria,  i.e.  experiment  altitude,  separation,  spatial  arrange¬ 
ment  of  the  satellites,  and  the  location  of  the  points  with  respect  to  the  sample 
points,  were  examined  for  geoid  heights  and  geoid  height  differences.  The 
results  for  geoid  height  differences  depend  mainly  on  the  level  of  correlation 
between  the  point  estimates.  For  small  distances,  e.  g.  150  km,  they  become 
as  low  as  0. 65  m  (altitude  200  km,  separation  250  km,  <7(  p)  =  ±0.03  mgal). 

For  point  geoid  heights  (and  the  same  situation  as  above)  about  ±0.80  m  are 
feasible.  In  addition,  the  information  content  in  the  three  frequency  ranges  from 
degree  13  to  120,  121  to  180,  and  181  to  infinity  was  investigated  separately.  The 
decrease  of  the  a  priori  variance  by  a  factor  of  12  in  the  frequency  range  from 
121  to  180  shows  that  sufficient  short-wavelength  information  -  down  to  half  wave¬ 
lengths  of  100  km  -  can  be  deduced  from  the  low-low  SST  experiment.  Throughout 
the  analysis  a  reference  field  up  to  degree  12  was  assumed  to  be  known  perfectly. 

Whereas  the  results  for  mean  gravity  anomalies  would  satisfy  in  resolution 
as  well  as  in  accuracy  the  defined  demands  of  geophysicists,  the  obtained  numbers 
for  geoid  heights  and  geoid  height  differences  are  not  sufficient  for  most  oceano¬ 
graphic  purposes.  On  the  other  hand,  it  is  worthwhile  to  remember  that  the  derived 
geoid  height  precision  of  ±0.80  m  is  very  well  comparable  with  the  excellent  result 
for  sea  surface  topography,  as  obtained  from  GEOS-3  altimetry. 


References 


Applications  of  a  Dedicated  Gravitational  Satellite  Mission,  National  Academy 
of  Sciences,  Washington,  D.C.,  1979. 

Barthel,  G. ,  T.  Halldorsson,  H.  Hufoagel,  D.  Meissner,  Ch.  Schmidt,  Laser 
Ranging  Instrumentation,  Final  Report,  Messerschmitt-Bolkow-Blohm, 

Space  Division  RX  12,  Munchen,  1978. 

Christodoulidis,  D.  C. ,  On  the  Realization  of  a  10  cm  Relative  Ocean  Geoid, 

Department  of  Geodetic  Science,  Report  No.  247,  The  Ohio  State  University, 
Columbus,  1976. 

Douglas,  B.  C.,  C.  C.  Goad,  F.  F.  Morrison,  Determination  of  the  Geopotential 
from  Satellite-to-Satellite  Tracking  Data,  NO  A  A  Techn.  Memo.  NOS  NGS 
24,  Rockville,  1980. 

Eddy,  W. ,  R.  Sutermeister,  Satellite-to-Satellite  Measurements,  Wolf  Research 
and  Development  Group,  Riverdale,  Maryland,  1975. 

Grafarend,  E. ,  The  Bruns  Transformation  and  a  Dual  Setup  of  Geodetic  Obser¬ 
vational  Equations,  preprint  NOAA  NOS  NGS  Report,  Rockville,  1979. 

Hajela,  D.  P. ,  Improved  Procedures  for  the  Recovery  of  5°  Mean  Gravity  Anomalies 
from  ATS -6 /GEOS -3  Satellite-to-Satellite  Range-Rate  Observations,  Depart¬ 
ment  of  Geodetic  Science,  Report  No.  276,  The  Ohio  State  University,  Colum¬ 
bus,  1978. 

Heiskanen,  W.  A.,  H.  Moritz,  Physical  Geodesy,  Freemann  and  Co. ,  San  Fran¬ 
cisco,  1967. 

Jekeli,  Ch. ,  Global  Accuracy  Estimates  of  Point  and  Mean  Undulation  Differences 
Obtained  from  Gravity  Disturbances,  Gravity  Anomalies  and  Potential  Coeffi¬ 
cients,  Department  of  Geodetic  Science,  Report  No.  288,  The  Ohio  State 
University,  Columbus,  1979. 

Kaula,  W.  M. ,  G.  Cherubini,  N.  Burkhard,  D.  D.  Jackson,  Applications  of 

Inversion  Theory  to  New  Satellite  Systems  for  Determination  of  the  Gravity 
Field,  Final  Report,  Department  of  Earth  &  Space  Sciences,  University  of 
California,  Los  Angeles,  1978. 

Krynski,  J. ,  Possibilities  of  Low-Low  Satellite  Tracking  for  Local  Geoid  Improve¬ 
ment,  Mitteilungen  der  geodatischen  Institute  der  TU  37,  Graz,  1978. 

Moritz,  H.,  Advanced  Least-Squares  Methods,  Department  of  Geodetic  Science, 
Report  No.  175,  The  Ohio  State  University,  Columbus,  1972. 


Muller,  P.  M. ,  W.  L.  Sjogren,  Mascons:  Lunar  Mass  Concentrations,  Science, 

161,  pp.  680-684,  1968. 

Rapp,  R.  H. ,  D.  P.  Hajela,  Accuracy  Estimates  of  l°x  1°  Mean  Anomaly 

Determinations  from  a  High-Low  SST  Mission,  Department  of  Geodetic  Science, 
Report  No.  295,  The  Ohio  State  University,  Columbus,  1979. 

Reed,  G.  B. ,  Application  of  Kinematical  Geodesy  for  Determining  the  Short  Wave 
Length  Components  of  the  Gravity  Field  by  Satellite  Gradiometry,  Depart¬ 
ment  of  Geodetic  Science,  Report  No.  201,  The  Ohio  State  University, 

Columbus,  1973. 

Rummel,  R. ,  Determination  of  Short- Wavelength  Components  of  the  Gravity  Field 
from  Satellite-to-Satellite  Tracking  or  Satellite  Gradiometry  -  An  Attempt 
to  an  Identification  of  Problem  Areas,  manuscripts  geodaetica,  4,  pp.  107- 
148,  1979. 

Rummel,  R. ,  D.  P.  Hajela,  R.  H.  Rapp,  Recovery  of  Mean  Gravity  Anomalies 
from  Satellite-to-Satellite  Range  Rate  Data  Using  Least  Squares  Collocation, 
Department  of  Geodetic  Science,  Report  No.  248,  The  Ohio  State  University, 
Columbus,  1976. 

Rummel,  R. ,  Ch.  Reigber,  K.  H.  Ilk,  The  Use  of  Satellite-to-Satellite  Tracking 
for  Gravity  Parameter  Recovery,  Proc.  European  Workshop  Space  Oceanog¬ 
raphy,  Navigation  and  Geodynamics,  ESA  SP 137,  pp.  153-161,  1978. 

Rummel,  R. ,  K.-P.  Schwarz,  M.  GerstI,  Least  Squares  Collocation  and  Regulari¬ 
zation,  Bull.  Geod.,  53,  pp.  343-361,  1979. 

Schwarz,  C.,  Gravity  Field  Refinement  by  Satellite-to-Satellite  Doppler  Tracking, 
Department  of  Geodetic  Science,  Report  No.  147,  The  Ohio  State  University, 
Columbus,  1970. 

Sunkel,  H. ,  A  Covariance  Approximation  Procedure,  Department  of  Geodetic  Science, 
Report  No.  286,  The  Ohio  State  University,  Columbus,  1979. 

Tscheming,  C.  C. ,  Covariance  Expressions  for  Second  and  Lower  Order  Deriva¬ 
tives  of  the  Anomalous  Potential,  Department  of  Geodetic  Science,  Report  No. 
225,  The  Ohio  State  University,  Columbus,  1976. 

Tscherning,  C.  C.,  R.  H.  Rapp,  Closed  Covariance  Expressions  for  Gravity 

Anomalies,  Geoid  Undulations  and  Deflections  of  the  Vertical  Implied  by  Anomaly 
Degree  Variance  Models,  Department  of  Geodetic  Science,  Report  No.  208, 

The  Ohio  State  University,  Columbus,  1974. 

Wolff,  M. ,  Direct  Measurements  of  the  Earth's  Gravitational  Potential  Using  a 
Satellite  Pair,  Journal  Geophys.  Res.,  14,  pp.  5295-5300,  1969. 


-42- 


Appendix 


The  degree  variance  model  underlying  all  further  derivations  is  the  one 
published  by  Tseheming  and  Rapp  (1974)  for  gravity  anomalies 


g^  +  3  A  (l-  1) 

(l- 2)(X  +  B) 


_  1*2 
s  c„ 


(Al) 


where  s  =  ( Ta^/r?2)  is  the  square  ratio  of  the  radius  of  an  adopted  Bjerhammar 
sphere  and  the  geocentric  distance  of  the  point  under  consideration,  A  =  425.  28 
mgafand  B  =  24.  Then  the  degree  variances  of  the  disturbing  potential,  and  of 
the  radial  and  horizontal  components  of  the  gradient  of  the  disturbing  potential 
become; 
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al<T)  = 

<  +  i_rsj» 

■  L ' 

(A2) 

»|(6r)  = 

S1*3  cl  .  and 

(A3) 

(A4) 

The  square  of  the  velocity  difference  of  two  satellites  is,  according  to  e.g.  Wolff 

(1969),  approximated  by  6xau  «  Tj  where  XeS  =  ("”)  is  the  mean 

velocity.  Thus  the  degree  variance  model  for  a  velocity  difference  for  instance 
in  radial  direction  becomes: 


°i  <!«**  i>  *  (*  *  (i^rT  °i « T|  <A5» 

with  rp  and  rQ  the  geocentric  radii  of  the  two  satellites.  Similarily,  the  degree 
variance  model  for  an  acceleration  difference,  again  in  radial  direction,  is: 


Equations  (A5)  and  (A6)  are  slightly  optimistic  degree  variance  models  for  range 
rate,  and  range  rate  change  signals.  They  are  optimistic  in  the  sense  that  one 
assumes  the  two  satellites  to  be  arranged  in  such  a  way  that  they  either  sense 
perfectly  the  radial  or  the  horizontal  component.  Finally,  the  degree  variance 
model  for  the  second  radial  derivative  is 


»!< T 


rr)  =  s**3 


(t-t-1)3  ( L  +  2)2 
r0„2  (l -If 
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(A7) 


A  comparison  of  signal  and  noise  requires  a  degree  or  degree-order  model 
for  the  noise,  too.  We  postulate  white  noise  behavior  and  approximate  its  covari¬ 
ance  function  by  a  first-order  Gauss-Markov  model,  e.g.  for  range  rates: 
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with  mo3  .. .  variance  of  the  random  observation  error, 
c  ...  inverse  correlation  length,  and 
0  ...  spherical  distance. 

A  harmonic  analysis  on  the  sphere  yields  for  a  coefficient  of  degree  l  and  order 
m  of  the  noise  process.- 

<|,(6)  -  t  ( A9i 

The  inverse  correlation  length  c  (units  of  arc*  is  derived  from  the  correlation 
length  in  time  units,  tc  ,  or  length  units,  sc ,  by: 


where  p  *  2^  (  rP3/GM)  is  the  orbital  period.  Equation  (A9)  holds  for  the  degree- 
order  variance  of  the  noise  of  any  type  of  observed  quantity  if  only  the  corresponding 
noise  variance  and  correlation  length  are  inserted. 


